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ABSTRACT 



Aims. Deviations from axial symmetry are necessary to maintain self-sustained MRI-turbulence. We define the parameters region 
where nonaxisymmetric MRI is excited and study dependence of the unstable modes structure and growth rates on the relevant 
parameters. 

Methods. We solve numerically the linear eigenvalue problem for global axisymmetric and nonaxisymmetric modes of standard-MRI 
in Keplerian disks with finite diffusion. 

Results. For small magnetic Prandtl number the microscopic viscosity completely drops out from the analysis so that the stability 
maps and the growth rates expressed in terms of the magnetic Reynolds number Rm and the Lundquist number S do not depend 
on the magnetic Prandtl number Pm. The minimum magnetic field for onset of nonaxisymmetric MRI grows with Rm. For given 
S all nonaxisymmetric modes disappear for sufficiently high Rm. This behavior is a consequence of the radial fine-structure of 
the nonaxisymmetric modes resulting from the winding effect of differential rotation. It is this fine-structure which presents severe 
resolution problems for the numerical simulation of MRI at large Rm. 

Conclusions. For weak supercritical magnetic fields only axisymmetric modes are unstable. Nonaxisymmetric modes need stronger 
fields and not too fast rotation. If Pm is small its real value does not play any role in MRI. 
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1. Introduction 

The leading concept for the origin of turbulence in accretion 
disks is the turbulence driving by the magnetorotational insta- 
bility (MRI). The instability can be excited by even a very weak 
magnetic field provided that rotation with outward decreasing 
angular velocity is present (see Balbus & Hawley [1998l ). 

MRI is expected to possess the remarkable property of being 
self-sustained, i.e. to support the destabilizing magnetic field via 
its own dynamo (Brandenburg et al. |1995l Hawley et al. 1996). 
The MRI-driven dynamo needs a finite but small initial field to 
start. Differential rotation alone cannot drive a dynamo (Elsasser 
1946). Deviations from axial symmetry are necessary for dy- 
namo and the deviations have to be produced by the MRI itself. 
The excitation of nonaxisymmetric modes of MRI is thus neces- 
sary for the self-sustained turbulence. 

The present paper focuses on the nonaxisymmetric MRI. A 
model of differentially rotating disk with finite diffusivities and 
axial background field is applied to perform linear analysis of 
global stability. As in the axisymmetric case, the nonaxisym- 
metric MRI can be found in a range between some minimum 
B m [ n and maximum Z? max values of the background field. The in- 
stability range depends on rotation rate (parameterized by the 
magnetic Reynolds number Rm). In contrast to the axisymmet- 
ric case, B m \ n for nonaxisymmetric modes does not approach a 
(small) constant value with increasing rotation but grows in lin- 
ear proportion to Rm. The larger Rm, the stronger background 
field is required to maintain the nonaxisymmetric instability. 



This behavior is increasingly difficult to follow numerically 
as Rm grows because of the winding effect of differential rota- 
tion. The shearing of nonaxisymmetric fields by differential ro- 
tation produces fine radial structure when the field is too weak 
to resist the winding. This shearing effect can simultaneously be 
the reason for the increase of B m in with Rm and for the high reso- 
lution required to resolve the nonaxisymmetric MRI numerically 
(Fromang & Papaloizou l2007l) . 



2. The model 



We use the model of Kitchatinov & Mazur d 19971 1 considering 
a rotating disk of constant thickness, 2H, threaded by a uniform 
axial magnetic field. The rotation axis is normal to the disk plane 
and the angular velocity, £2, varies only with the distance s to the 
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This profile describes almost uniform rotation with the veloc- 
ity £2o up to the distance so and Kepler rotation for larger dis- 
tances. The aspect ratio is s^/H = 5. We assume incompress- 
ibility, div U = 0. The motion equation is curled to exclude the 
pressure. This results in the equation for vorticity, W = V x U, 



dW 

— = Vx(UxW + Jx B/p) + vAW, 
dt 



(2) 
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where B is magnetic field, J = n Q l V x B is the current density, 
and v is the viscosity. The magnetic field follows the induction 
equation 



dB 

— = V X (U X B - nV X B) . 

at 



(3) 



The equations are linearized about the background state of ro- 
tational motion (fTJ and uniform axial magnetic field Bq = z.Bq, 
where z is unit vector along the rotation axis. Pseudovacuum 
conditions for the magnetic field, z x B' = 0, are applied on the 
disk surfaces at z = ±H, where prime signifies small disturbance. 
The conditions for the flow assume impenetrable and stress-free 
boundaries. To treat large radial distances we change to a new 
variable, 



y = 



1 + s/so' 



0<y < 1. 



(4) 



and use uniform numerical grid in y. In effect, we did not com- 
pute up to y — 1 but imposed a side boundary at y — 0.99 with 
the condition that all disturbances vanish at this boundary. 

The magnetic and velocity disturbances are expressed in 
terms of scalar potentials, e.g. 



B' = z x Vfi + V x (z x VA) , 



(5) 



so that the disturbances are automatically divergence-free. The 
Fourier expansions in z and in the azimuthal coordinate of the 
cylinder coordinate system, 

B = e rt J] e im ^ (Bfjs) cos - 1 /2)z) + Bf m (s) sminlz)) , (6) 



are applied. Summation in Eq. (|6]l runs over / = 1, 2, 3, ... and 
m = 0, 1,2, .... The mathematical treatment of the velocity and 
vorticity disturbances are quite similar. 

The equation system for the disturbances split into a set of in- 
dependent equations for different / and m. Also the terms marked 
by the upper indeces S and Am© are not mixed by the equa- 
tions. The indices signify the magnetic modes symmetric and 
antisymmetric relative to the midplane of the disk. We shall use 
the notation Sm and Am for the symmetric and antisymmetric 
modes where m is the azimuthal wave number. Note that Sm and 
Am represent families of modes that can be further distinguished 
by the vertical wave number /. For fixed m, I and the symme- 
try type we have an eigenvalue problem for ordinary differential 
equations in the variable y. The problem is solved numerically. 

The problem has three governing parameters, i.e. the mag- 
netic Reynolds number (Rm), the Lundquist number (S), 



Rm = 



Q {) H 2 



s = 



VaH 



(7) 



where Va = Bo/ yfjjp is the Alfven velocity, and the magnetic 
Prandtl number 



Pm = -. 



(8) 



The Reynolds number Re = Rm/Pm and the Hartmann number 
Ha = S/ VPm were other possible choices. We shall see, however, 
that the parameters set (0 is most convenient for MRI analysis 
for small Pm. 



3. Marginal stability from local analysis 

For disturbances of small spatial scale compared to local radius, 
the differential rotation can be approximated by a plane shear 
flow which leads to the local approximation (Hawley & Balbus 
|19911 l. For the simplest case of plane-wave disturbances with k\\z 
one finds the dispersion relation 

(cr + J]k 2 f ((cr + vk 2 f + 2 (2 - q) £2 2 j + 

+Q A (Q A - 2qQ 2 + 2 (cr + vk 2 ) (cr + rjk 2 )) = 0, (9) 

where cr is the eigenvalue, £2 A = kV\ is the magnetic fre- 
quency and q — -d log O/d log s is the local shear. An instability 
for the considered case can come only as a change of stability 
(Chandrasekhar [T96 1 l l . Therefore, we can put cr — in Eq. (O to 
get the neutral stability equation defining the boundary between 
stability and instability (Kitchatinov & Riidiger 2004): 



Rm 2 = 



(Pm + S 2 ) 2 
2{qS 2 -2 + q)' 



(10) 



The dimensionless quantities © are redefined in terms of the 
wave number (H — > k^ 1 ). 

Equation ( fTOb shows that the instability requires sufficiently 
large Rm exceeding 



Pm^(l + ^ 
q \ qPm 



(11) 



which corresponds to S = (Pm + 2(2 - q)/q) 1 ^ 2 . For Rm » 
Rm m j n , the instability only exists for S between a lower and an 
upper limit, i.e. J(2 - q)/q < S < -\j2~q Rm. For small Pm the 
Eqs. ([Tol l and ( fTTT i loose their dependence on Pm so that the vis- 
cosity completely drops off. The instability in this limit is fully 
controlled by Rm and S both defined in terms of the magnetic 
diffusivity. The computations of the most unstable axisymmet- 
ric modes of MRI for our global model confirm these results. 
FigureQ] shows that the entire stability map in the Rm-S plane 
does not depend on Pm for small Pm (Riidiger & Kitchatinov 
120051 1 . Also the growth rates of the linear instability for small 
Pm are fully controlled by Rm and S . 



4. Nonaxisymmetric modes 

The magnetic Prandtl number is very small for cool protostel- 
lar and protoplanetary disks, it is also small for liquid metals 
of laboratory experiments but it is very large for galaxies (cf. 
Brandenburg & Subramanian 2005). Properties of MRI can be 
expected to vary strongly between the two cases of small and 
large Pm (see Lesur & Longaretti 2005; Fromang et al. 120071 1. 
The linear theory, however, shows that for Pm < 1 the viscos- 
ity does not play any role. This is in contradiction with Lesur & 
Longaretti (2007). 



4.1. Small Pm 

The results for the nonaxisymmetric modes confirm that MRI for 
small Pm is not sensitive to viscosity. The neutral stability lines 
for m — 1 given in Fig.|2]approach one and the same contour as 
Pm decreases. The growth rates show the same tendency (Fig.[3j. 
Note that one cannot distinguish the cases of Pm = 0.01 and Pm 
= 0.001 by their growth rates or stability maps. 
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Fig. 1. Neutral stability lines for most easily excited SO modes 
and various Pm< 1 . The lines do not depend on Pm provided that 
Pm is small. 
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Fig. 2. Neutral stability lines for nonaxisymmetric S 1 modes of 
the lowest vertical wave number (/ = 1) and various Pm. The 
lines are marked by the corresponding magnetic Prandtl num- 
bers. 



This simple scaling should be important for numerical sim- 
ulations. The magnetic Prandtl numbers in astrophysical bodies 
can be so small that they are not accessible for simulations. We 
have shown with our linear theory that the basic instability for 
moderately small magnetic Prandtl number (Pm ^0.1) closely 
reproduce the results for indefinitely small Pm (provided that the 
results are expressed in terms of Rm and S or other parameters 
not including the viscosity). The scaling means that MRI at small 
Pm does not develop so fine structure that the viscosity becomes 
important. 

The strong-field limit of the instability domain in Fig. ^be- 
haves like Rm ^ 2 S, i.e. the rotation is slightly super Alfvenic. 
The new feature in Fig. [2] is that the minimum field, B m ; n , pro- 
ducing the instability also grows with Rm. This positive slope is 
in contrast to the axisymmetric modes of Fig. Q] for which B m ; n 
approaches a small constant value as Rm grows. For given field 
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Fig. 3. Growth rates of nonaxisymmetric S 1 modes of the lowest 
vertical wave number (1=1) for various Pm. The plot shows 
the growth rates in the range of magnetic fields producing the 
instability forRm = 150. 

strength, the nonaxisymmetric modes are switched off by too 
fast rotation. The relation Rm st 10 S provides an estimation 
for the upper limit on rotation rate for the S 1 modes excitation. 
For Rm > 10 S the SI modes with / = 1 no longer exist. With 
a plasma-/? with ft =s Rm 2 /S 2 one finds f3 S 100 as the insta- 
bility condition for S 1 modes with I = 1 . Hence, the magnetic 
field must be strong enough. We shall demonstrate below that for 
Rm > 10 S the realization of the mode stability requires rather 
high numerical resolution in radial direction. The nonaxisym- 
metric modes are necessary for self-sustained MRJ-turbulence. 
A self-sustained turbulence in high Rm-regime may thus not be 
possible. Another possibility is the nonaxisymmetric instability 
of an imposed azimuthal magnetic field ("AMRI", Riidiger et al. 
120071 Simon & Hawley |2009l l. 

The increase of B m i n with Rm for nonaxisymmetric MRI 
is a consequence of the winding effect of differential rotation. 
The pitch-angle of unstable disturbances near fi m j n is small 
(Kitchatinov & Riidiger 2004), i.e. the winding is strong. The 
differential rotation converts azimuthal inhomogeneity of the 
nonaxisymmetric modes into a fine radial structure which is fi- 
nally destroyed by diffusion. 

The increase of B m [ n with Rm also appears for large Pm. In 
this case it is reasonable to use Re and S/Pm to parameterize the 
rotation and the background field. The lines show little depen- 
dence on large Pm when plotted in the plane of these parameters 
which now do not depend on the magnetic diffusion 77. 

4.2. Overtones 

Another new feature of the nonaxisymmetric instability is that 
modes with higher vertical wave number (I > 1) are preferred 
for some parameter domains. For axial symmetry, the region of 
parameters where the SO mode with / = 1 is unstable includes 
the instability regions of all other modes. 

Figure|4] shows the neutral stability lines for several non- 
axisymmetric modes together with the line for the most unsta- 
ble axisymmetric SO mode. The lines for the nonaxisymmetric 
modes intersect so that the modes with finer vertical structure 
are preferred on the weak-field side of the stability map. This is 
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Fig. 4. Common stability map of SO mode (/ = 1) and several 
nonaxisymmetric modes of different structure in z-direction. S 1 
mode is preferred on the strong-field side of the plot. On the 
weak-field side, the neutral stability lines of nonaxisymmetric 
modes intersect. Pm = 1 . 




Fig. 5. Vectorplots of the unstable SI modes (I — 1) for the mid- 
plane of the disk after computations for Rm =100 and Pm = 1. 
The left panel is for the field slightly exceeding B m \ a and the 
right panel - slightly below B max . The dashed circle shows the 
corotation radius s — sq. 



again related to the winding effect of differential rotation. The 
modes with finer vertical structure have larger Lorentz force to 
resist the winding of weak fields. Note that only the SI mode 
with / = 1 can compete with SO mode. There is a narrow region 
on the strong-field side of Fig.|4]where the mode SO is stable but 
SI not. 

The positive slope on the upper Rm(S) curve of the stability 
map also exists if higher / modes are included. The existence of 
nonaxisymmetric magnetic instability is necessary for any form 
of MRI-dynamo. The preference of nonaxisymmetry for strong 
fields is promising for the dynamo and the self-sustained tur- 
bulence. The field strength should increase sufficiently fast with 
Rm, or the plasma-/? of shearing box simulations (Fromang et al. 
2007) should not be too high, in order to probe this possibility. 
Otherwise, for too large B, i.e. for too weak fields the nonax- 
isymmetric modes do not appear. 
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Fig. 6. Growth rates of SI modes for Re = 150, Pm = 1 (left), 
and Pm = 10 (right) computed with different resolutions. The 
dotted, dashed, full, and dashed-dotted lines were obtained with 
100, 150, 200, and 300 grid points in radius, respectively. 

5. The resolution problem 

The upper branches of the neutral stability lines for the nonax- 
isymmetric modes in Fig. [4] are always terminated because of 
numerical resolution problem. The needed resolution to follow 
the lines for higher Rm rapidly increases. Again the shearing by 
differential rotation is the reason. Figure|5] gives vector plots of 
unstable nonaxisymmetric magnetic modes for weak (close to 
B m \ n ) and strong (close to B max ) background fields. Obviously, 
the disturbances in the strong-field case resist the shearing by 
the differential rotation. On the other hand, the shearing of weak 
fields produces tightly wound spirals with increasing require- 
ments for the numerical resolution. 

In Fig. [6] the MRI growth rates calculated for the SI mode 
computed with different numbers of radial grid points are shown. 
There is no resolution problem at all for strong fields. For weak 
fields, however, a too low resolution produces an unreal insta- 
bility. This numerical artifact can be suppressed by increasing 
resolution. For fixed Reynolds number it is even harder to do it 
for large Pm. In the given examples always more than 200 grid- 
points are necessary to overcome the numerical artefacts. 

The resolution problem has also been met in nonlinear shear- 
ing box simulations of MRI (Fromang & Papaloizou 2007J. If 
our interpretation of the problem as a result of rotational shear- 
ing is correct then an increase of resolution only in one radial 
direction is necessary to resolve the problem. Also their value 
of the plasma-/? of order 400 leading to the relation Rm =s 20 S 
appears to imply too weak magnetic fields outside of the region 
of nonaxisymmetric MRI. 

We conclude that a self-sustained MRI-turbulence may only 
be found with sufficiently strong initial field. The minimum field 
for the turbulence exceeds by at least one order of magnitude the 
minimum external field required for axisymmetric MRI and it 
increases with magnetic Reynolds number. 
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